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Abstract We review recent progress on the lattice simulations of the chiral mag- 
netic effect. There are two different approaches to analyze the chiral magnetic effect 
on the lattice. In one approach, the charge density distribution or the current fluc- 
tuation is measured under a topological background of the gluon field. In the other 
approach, the topological effect is mimicked by the chiral chemical potential, and 
the induced current is directly measured. Both approaches are now developing to- 



p ; 

ward the exact analysis of the chiral magnetic effect 



1 Introduction 

In the strong interaction, the gauge field forms nontrivial topology. The existence of 
the topology has been theoretically established, while its observation is difficult in 
experiments. The chiral magnetic effect is a possible candidate to detect the topolog- 
ical structure in heavy-ion collisions HI. The chiral magnetic effect is the generation 
of an electric current in a strong magnetic field. 

The essence of the chiral magnetic effect is the imbalance of the chirality, i.e., 
the number difference between the right-handed and left-handed quarks. The mag- 
netic field induces the electric currents of the right-handed and left-handed quarks 
in opposite directions. If the chirality is imbalanced, a nonzero net electric current 
is induced. In a local domain of the QCD vacuum, the chiral imbalance is generated 
by the topological fluctuation and the axial anomaly. In the global QCD vacuum, the 
chirality is balanced as a whole. The strong theta parameter is experimentally zero, 
0=0, although its reason is unknown. This is the strong CP problem. The chiral 
magnetic effect is regarded as the local violation of the CP symmetry. 
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Fig. 1 A cartoon of how to observe the chiral magnetic effect on the lattice. Left: a topological 
charge Q of the gauge field induces a nonuniform current density distribution. Right: a chiral 
chemical potential /is induces a uniform electric current. 



Experimental facilities tried to measure the chiral magnetic effect through charged- 
particle correlations l2][3l. However, the interpretation of the experimental data is 
not yet conclusive. On the theoretical side, the chiral magnetic effect has been stud- 
ied in various frameworks, e.g., phenomenological models, the gauge-gravity dual- 
ity, etc. The chiral magnetic effect has been also studied in the lattice simulations. 
The lattice simulation is a powerful framework to solve QCD nonperturbatively on 
computers. By means of the lattice simulation, we can study the chiral magnetic 
effect from first principles in QCD. 

There are two approaches to analyze the chiral magnetic effect in the lattice simu- 
lation. In other words, there are two different ways to generate the chiral imbalance: 

1. topological charge il|5l|6l HE m 

2. chiral chemical potential ifTOl fTTI [Till 

These concepts are schematically depicted in Fig. [T] In the first case, a topologi- 
cal charge of the background gauge field generates the chiral imbalance, which is 
spatially nonuniform. When an external magnetic field is applied, a current density 
distribution appears around the topological object. In the second case, a chiral chem- 
ical potential generates the chirally imbalanced matter, which is spatially uniform. 
A uniform electric current is induced by the external magnetic field. 

In this chapter, we overview the theoretical background and the current status of 
the lattice studies of the chiral magnetic effect. Here we focus only on the lattice 
aspect of the chiral magnetic effect. For the theoretical and phenomenological as- 
pects, see the corresponding chapters. We use the Euclidean metric and the lattice 
unit in the following sections. 
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2 Basics of the lattice simulation 

The basic formalism of the lattice simulation has been well established. For the 
details, see the textbooks lfT3lfT4lll5l[T6l . The formalism is based on the Euclidean 
QCD partition function 

Z)t/detZ)[t/]e"'^Gl^l . (1) 



The space-time is discretized as a hypercubic lattice. The gluon field is written as 
the SU(3) Unk variable 

f/^(x)=expfef"A«(x)) . (2) 

The functional integral is numerically evaluated by the Monte Carlo simulation. We 
generate gauge configurations, which are sets of the link variable, and then calculate 
the expectation value of an operator as 

io[u]) = -^j:m- (3) 

JVconf {uj 

The gauge configurations are generated to satisfy the probability weight P = 
detD[f/]e^'^G['^]. The simulation including the fermion determinant is called the dy- 
namical QCD simulation or the full QCD simulation. The quenched approxima- 
tion is often used to reduce the computational cost. In the quenched approximation, 
the fermion determinant is ignored and the probability weight is P = e^^oW]^ jyiq 
quenched gauge configurations are independent of the fermion action. 

The probability weight must be positive real, otherwise it cannot be interpreted 
as the probability weight. In QCD, the fermion action becomes complex at a finite 
quark chemical potential. The definition of the probability weight must be modified, 
e.g., by the reweighting method ifTTl . Even after the modification, the Monte Carlo 
simulation severely suffers from strong sign fluctuation. This is known as the sign 
problem. The sign problem at the quark chemical potential is an important unsolved 
problem in the lattice simulation [18|. As shown later, a chiral chemical potential 
does not cause the sign problem. This is similar to two-color QCD |fT9l [2011211 and 
an isospin chemical potential 11211 l22l |23l l24l |25 l . 

The basic observable of the chiral magnetic effect is the local vector current 
density 

jii{x) = W{x)7ti¥{x) ■ (4) 



The fourth (zeroth) component corresponds to the local charge density. For calcu- 
lating the local vector current density, we consider the Dirac eigenvalue problem 
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D[U](l)k{x)^{ih + m)Ux) , (5) 

and use the identity 

{wi^)7M^)) - . . (6) 

Thus, its expectation value is obtained by inverting or diagonalizing the Dirac op- 
erator D[U]. In the case of diagonalizing, we can calculate the local vector current 
density of each Dirac eigenmode. 

Chiral symmetry is a nontrivial problem on the lattice due to the Nielsen- 
Ninomiya no-go theorem 1261 |27l . Most lattice fermions more or less break chiral 
symmetry. The lattice fermion with exact chiral symmetry has been known, although 
its computational cost is rather large in the dynamical simulation. We should select 
an appropriate lattice fermion, corresponding to the purpose of the simulation. For 
the details of the lattice fermions and chiral symmetry, see the reviews ||28ll29l[30l . 

For a magnetic field, the QED gauge field is also introduced. When the magnetic 
field is external, i.e., not dynamical, the QED field strength term does not exist in the 
action. To couple the fermions to the magnetic field, the Dirac operator is replaced 
as 



D[U]^D[uU] (7) 

with the U(l) link variable 

Ufi{x) = exp(/g'A^(x)) . (8) 

We can apply a homogeneous magnetic field in a finite-volume box with periodic 
boundary conditions. For the homogeneous magnetic field in the z-direction, the 
U(l) hnk variables are set as 

ui{x) ^ exp{-iqBNsy) atx^Ns (9) 
U2{x) = exp{iqBx) (10) 
u^{x) ^ I for other components (11) 

in the lattice volume A'^"' x Nt ISTl . In this setup, the magnetic field is quantized as 

2n 

= X (integer) . (12) 



This integer is the input parameter which controls the strength of the magnetic field 
in the simulation. 
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3 Lattice simulation with a topological background 

The gauge configuration possesses a topological charge. The topological charge of 
the gauge configuration is given as 

Euclidean topological objects, such as the instanton, can be reproduced on the lattice 
when the gauge configuration is smooth enough [32]. 

The fermion feels the background topology of the gauge configuration through 
the zero mode. The zero mode is defined as the eigenmode 0^ which has the zero 
eigenvalue iXk = in Eq. The topological charge and the zero mode are related 
through the Atiyah-Singer index theorem 

NR-NL^NfQ, (14) 

where A^^ and A^^ are the numbers of the right-handed zero modes and the left- 
handed zero modes, respectively i33|. The fermion zero mode is essential to gen- 
erate the chiral imbalance in the topological background. We must use the lattice 
fermion which is sensitive to the zero mode and satisfies the index theorem, e.g., the 
overlap fermion. 

Naively, it is impossible to measure the local vector current density ^ in the 
topological background. The reason is as follows. In the QCD vacuum, the positive 
and negative topological charges appear with the same probability. In the simulation, 
the numbers of the gauge configurations with the positive and negative topological 
charges are the same, namely 

(e) = o. (15) 

The positive and negative topological charges induce the vector current in opposite 
directions. When all the Monte Carlo samples are averaged in all the topological 
sectors, the net vector current is zero. To measure the vector current, we must fix the 
topological sector by the lattice action which suppresses topology changing tran- 
sitions 1341 . Although the fixed-topology simulation cannot reproduce the 0=0 
vacuum, we can obtain a finite expectation value of the vector current. 

The fixed-topology analysis has been done in the (2H-l)-flavor dynamical QCD 
simulation with the domain-wall fermion HHH. This simulation includes not only 
the external magnetic field but also the dynamical QED effect. The domain-wall 
fermion does not have the exact zero mode due to small explicit chiral symmetry 
breaking, but has the "near" zero mode which becomes the exact zero mode in an 
ideal limit. In Fig. 121 we show the charge density distribution of one near zero mode 
in one typical gauge configuration |[8]. The charge density of the A:-th eigenmode is 
defined as 
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Pkix) = , — ■ (16) 

lAk + m 

The simulation was performed above the critical temperature. The charge density 
distribution at Z? 7^ differs from that at Z? = 0. This result suggests that some relation 
exists between the charge density and the magnetic field. However, the exact relation 
is not clear in this simulation. We need further investigation for evidence of the 
chiral magnetic effect. 




Fig. 2 The charge density distribution p{x) in the (2+l)-flavor dynamical QCD+QED simulation 
at a^qB = (left) and 0.0490874 (right) [S]. The magnetic field is applied along the z-axis. The 
temperature is above the critical temperature. 



Except for the fixed-topology simulation, the vector current itself is zero because 
of the CP oddness. In this case, a numerical observable is the CP-even quantity 
which reflects the topological fluctuation 

~ (200 MeV)"* . (17) 

For instance, the fluctuation of the vector current is CP-even. This situation is sim- 
ilar to the experimental observation. An experimental observable must be CP-even, 
although the chiral magnetic effect is a CP-odd process. We have to extract the CP- 
odd information from the CP-even particle correlation. This kind of analysis is not 
easy because the fluctuation can be easily induced by other irrelevant effects. The 
irrelevant contributions must be subtracted correctly. 

The fluctuation (y^) of the vector current was calculated in the quenched SU(2) 
simulation at zero temperature lH, in the quenched SU(2) simulation at finite tem- 
perature 13] m, and in the quenched SU(3) simulation fT\. The overlap Dirac oper- 
ator was adopted in these simulations, although the zero modes were ignored. The 
vector currents are zero in all the directions because the topological sector is not 
fixed, but the current fluctuation is nonzero. In Fig.|3] we show the current fluctua- 
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tion in the quenched SU(2) simulation below and above the critical temperature Q. 
The ultraviolet part of the fluctuation is subtracted to obtain a clear signal as 

Ul)lR = ^T.{JI{x))b.T -^Y.{jlix))B=Q.T=0 , (18) 
site site 

where the index jj. is not summed over At zero temperature T = 0, all the fluctua- 
tions grow at stronger magnetic field. In particular, the longitudinal fluctuation {j^) 
grows faster than transverse fluctuations (y'j ) — Above the critical temperature 
T > Tc, the longitudinal fluctuation is insensitive and the transverse fluctuations de- 
crease at stronger magnetic field. As a consequence, the ratio of the longitudinal 
fluctuation to the transverse fluctuation is enhanced by the magnetic field in both 
cases. 
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Fig. 3 Tlie current fluctuation (i^)iR in the quenched SU(2) simulation at T = (left) and T 
1.1271. (right) El. The magnetic field is applied in the = 3 direction. 



As shown above, the magnetic field affects the charge density distribution and the 
current fluctuation. Note however that we must carefully check whether its origin is 
actually the chiral magnetic effect. In general, a strong magnetic field can induce a 
strong current fluctuation in the longitudinal direction, even if there is no topological 
object. This complication is the same as that in experiments. For identifying the 
chiral magnetic effect, we must distinguish a small topological contribution from 
other large contaminations in a high-precision simulation. 



4 Lattice simulation with a chiral chemical potential 

Another possible source of the chiral imbalance is a chiral chemical potential. The 
chiral chemical potential jLi5 is defined as 

D{jXs) = (<5m + igt"A"^ {x)) + m + Ai57475 (19) 
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in the continuum space 1351 . The chiral chemical potential directly couples to the 
chiral charge 



Ns =Nr-Nl ^- J d^x{\f/{x)'Y4'Y5\i/{x)) . 



(20) 



By using the chiral chemical potential, we can generate a chirally imbalanced QCD 
matter in equilibrium. The chiral chemical potential is the external parameter which 
tunes the chiral charge instead of the topological charge. The chiral chemical po- 
tential does not exist in the original QCD action because the chiral charge is not a 
conserved quantity. It is not a "chemical potential" in the exact sense. 

Because the topological charge is not necessary in this approach, the sensitivity 
to the zero mode is not important for the choice of the fermion action. For example, 
the lattice Dirac operator of the Wilson fermion is 

-D^i^s) = l^KY,[il- 7i)Ti+ + (1 + Y,)T,J\ 



m 



-ic[(l ~ 74e^'^')T4+ + (1 + y4e-^>^')T4-] , (21) 



with 



1 



(22) 
2m + 8 

[T^+],,y = U^{x)5,+ii,y (23) 
[T^-ky = U;{y)5,^f,,y. (24) 

The chiral chemical potential is introduced as the exponential matrix factor 

e±M5y5 ^ coshjLi5 ± 75 sinh^s , (25) 



which is the straightforward analogy to a quark chemical potential 11361 . The Wilson- 
Dirac operator (121]) reproduces the continuum form ( fT9] l in the continuum limit. 

A notable feature of the chiral chemical potential is that it does not cause the 
sign problem unlike the quark chemical potential. The Wilson-Dirac operator i2T[ 
is "75-Hermitian", 

Y5D{fi5) = [Y5D{ii5)]^ or 750(^^5)75 = £>'(Ai5) • (26) 
In the two-flavor case, the fermion determinant is positive real, 

det f ^^^^^ nn, ^ ) = det/)(M5)det75Z)(M5)75 = \deW{^5)\^ > . (27) 



V ^(A's) 

Therefore there is no sign problem. We can exactly simulate a kind of finite density 
QCD matter by the chiral chemical potential. 
In Fig.m we show the chiral charge density 
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«5 = ^ = -^L(V^Wr475rW) (28) 

' ' site 

of the Wilson fermion in the two-flavor dynamical QCD simulation |[T2| . The lattice 
spacing is a ~ 0.13 fm. The physical temperature is T ^ 400 MeV, which is above 
the critical temperature. The chiral charge density is finite at a finite chiral chemical 
potential. This means that the uniform chirally imbalanced matter is realized on 
the lattice. The total chiral charge in this lattice volume is A^5 = n^V ~ (9(10^). 
This number is much larger than a typical number of the topological charge. The 
typical number of the topological charge is (9(10) at most in the conventional lattice 
simulation. Owing to the large chiral imbalance, the analysis of the chiral magnetic 
effect becomes easy. 

0.4 



0.3 
^ 0.2 

m 
C 

0.1 


0.1 0.2 0.3 0.4 0.5 

Fig. 4 Tlie cliiral cliarge density nj in tlie two-flavor dynamical QCD simulation 1121 . The lattice 
spacing is a ~ 0. 13 fm and the temperature is T ~ 400 MeV. 



For the analysis of the chiral magnetic effect, the local vector current density (|4]i 
was measured. The vector current is induced only in the longitudinal direction of 
the magnetic field. The transverse components are exactly zero, (ji) = (72) = 0. In 
Fig.|5] the induced current 

site 

is plotted as a function of the magnetic field and of the chiral chemical potential. 
This induced current is direct evidence of the chiral magnetic effect. The induced 
current is a linearly increasing function in both cases. Therefore, the functional form 
is 




J = NiofCjXsqB 



(30) 
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Because all the fermions have the same charge in this simulation, the prefactor is 
Niof = 3 (color) X 2 (flavor) = 6. The overall coefficient C characterizes the strength 
of the induced current. This functional form is consistent with the analytical for- 
mula, 

-/ = A^dof TrrMs^fi , (31) 

which was derived from the Dirac equation coupled with the background magnetic 
field 135]. Note that Eq. ([3B is different from Ref. ||35l by q due to the definition 
of the electric current, i.e., /em = iJ- If there are several fermions with different 
charges, the total electric current is Jem = Hi^iJi = HiljCflsB. 




Fig. 5 The vector current density / in the two-flavor dynamical QCD simulation 1121 . The data are 
plotted as a function of the magnetic field B (left) and of the chiral chemical potential (right). 
The lattice spacing is a ~ 0. 13 fm and the temperature is 7 ~ 400 MeV. 



The overall coefficient is C = 0.013 ±0.001 in this lattice simulation and C = 
l/(2;r^) ~ 0.05 in the analytical formula. The induced current seems somehow 
smaller than the analytical formula. However, these overall coefficients should 
not be compared naively. For the quantitative argument, it is necessary to esti- 
mate several systematic effects in the lattice simulation. One important effect is 
the renormalization of the local vector current. The local vector current (|4|i is not 
renormalization-group invariant on the lattice 1371 . This property is different from 
that in the continuum theory. The local vector current is renormalization-group in- 
variant in the continuum theory because of the Ward identity. We must take the 
continuum limit to compare the induced currents on the lattice and in the contin- 
uum. By taking the continuum limit, we can also remove other lattice discretization 
artifacts. For example, the Wilson fermion expUcitly breaks the chiral symmetry due 
to the discretization artifact. 

The systematic errors were partly estimated in Ref. lT2l . By varying simula- 
tion parameters, the dependences of the overall coefficient were examined in the 
quenched simulation. Although the dynamical QCD simulation is necessary for the 
quantitative argument, the quenched simulation is useful to understand which sys- 



Chiral Magnetic Effect on the Lattice 



11 



tematic effect is important. Actually, the quenched results were quahtatively similar 
to the dynamical QCD results. It turned out that the induced current is insensitive 
to the temperature, the quark mass, and the spatial volume. However, the overall 
coefficient strongly depends on the lattice spacing. The overall coefficient increases 
near the continuum limit. This systematic analysis indicates that the continuum ex- 
trapolation is necessary for the quantitative argument. 

Another important effect is chiral symmetry. It is difficult to discuss chiral sym- 
metry using the naive Wilson fermion. The Wilson fermion explicitly breaks chiral 
symmetry at a finite lattice spacing, while the explicit breaking vanishes in the con- 
tinuum limit. One possible origin of the strong lattice spacing dependence might 
be this artificial chiral symmetry breaking. We should investigate the role of chi- 
ral symmetry in the chiral magnetic effect by performing the same analysis with a 
chiral lattice fermion, such as the domain-wall fermion or the overlap fermion. 



5 Conclusion 

In this review, we have overviewed the lattice studies of the chiral magnetic effect. 
The vector current and its fluctuation were measured in the chiral imbalance, which 
is generated by the topological charge or the chiral chemical potential. We should 
develop these pioneering works in future. In the future works, it is important to 
respect the essential pieces of the chiral magnetic effect, in particular, the fermion 
zero mode and chiral symmetry. 

We see that the chiral magnetic effect is an observable phenomenon on the lattice. 
The lattice simulation is a hopeful approach to study the chiral magnetic effect in 
"numerical" experiments. 
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